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Abstract 

We investigate the behavior of a pair of heavy fermions, denoted by Q and Q, 
in a hot/dense medium. Although we have in mind the situation where Q and Q 
denote heavy quarks, our treatment will be limited to simplified models, which bear 
only some general similarities with QCD. We study in particular the limiting case 
where the mass of the heavy fermions is infinite. Then a number of results can be 
derived exactly: a Schrodinger equation can be established for the correlator of the 
heavy quarks; the interaction effects exponentiate, leading to a simple instantaneous 
effective potential for this Schrodinger equation. We consider simple models for the 
medium in which the QQ pair propagates. In the case where the medium is a plasma 
of photons and light charged fermions, an imaginary part develops in this effective 
potential. We discuss the physical interpretation of this imaginary part in terms of 
the collisions between the heavy particles and the light fermions of the medium; the 
same collisions also determine the damping rate of the heavy fermions. Finally we 
study the connection between the real-time propagator of the heavy fermion pair 
and its Euclidean counterpart, and show that the real part of the potential entering 
the Schrodinger equation for the real-time propagator is the free energy calculated 
in the imaginary-time formalism. 
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1 Introduction 



The propagation of heavy quark bound states, such as the J/^f, in a hot 
and dense medium has been the subject of many studies, since the original 
proposal by Matsui and Satz [1] suggesting that the suppressed production of 
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such bound states in ultra-relativistic heavy-ion collisions could be used as a 
signal of deconfinement. The initial suggestion [1] assumed that the dominant 
effect of the interaction with the medium is the screening of the QQ interaction 
potential, leading to a disappearance of the bound states when the screening 
is sufficiently strong. The screened potential can be obtained from equilibrium 
calculations of the force between two infinitely heavy color charges, a quantity 
which can, in principle, be extracted from the lattice data for the correlator 
of two Polyakov lines. Such a correlator is in fact related to the change of free 
energy occurring when a heavy QQ pair is added into a hot QCD medium 
[2], and the relation between this free-energy and the effective potential to be 
inserted in a Scrodinger equation (or in a T-matrix calculation, like in [3,4]) 
remains an open question. Thus, the effective potential has been identified 
alternatively with the QQ free energy [5-8], the internal energy [9,10] or even 
with a linear combination of both [11-13], leading to different dissociation 
temperatures for the J/^ff, depending on the choice made. 

In another line of investigation one considers the spectral functions of a cc 
and bb mesons in a hot environment [14-18]. Such objects contain in principle 
the full physical information on the properties of the mesonic excitations in a 
medium, and are also directly related to observable quantities, like the dilepton 
production rate in the vector channel. Unfortunately spectral functions are not 
directly measured in a lattice simulation, but have to be reconstructed from 
Euclidean correlators, and this reconstruction is plagued with large uncer- 
tainties. To go around such difficulties, one may attempt to compare directly 
the spectral functions and the corresponding imaginary-time correlators that 
are obtained from lattice calculations to those that can be calculated from a 
Schrodinger equation with a screened potential. Such studies were carried out 
recently [4,13,19-22]. 

More recently, still another approach was developed, in which the focus is 
on the real-time propagator of a heavy-quark pair. In Refs. [23-27], it was 
argued that such a propagator obeys a Schrodinger equation with an effective 
potential which was derived first perturbatively, and then estimated in classical 
lattice gauge simulations. An interesting observation made in [23] is that this 
effective potential may contain an imaginary part that may strongly affect the 
calculated spectral function. 

In this paper we shall discuss a number of issues triggered by the latter stud- 
ies. We shall examine under which conditions, the in-medium QQ-propagator 
obeys a Schrodinger equation. To do so, we shall consider the general equation 
of motion satisfied by the heavy fermion pair correlator. This is part of an infi- 
nite hierarchy of coupled equations for the n-point functions. However, in the 
limit where the mass of the heavy fermion is infinite, the equation for the QQ 
correlator decouples and reduces to a closed, Schrodinger equation. We shall 
also establish a link beween the effective potential entering this Schrodinger 
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equation and the heavy-quark free-energy that can be extracted from lattice 
data; in other words, we shall investigate the relations between the real-time 
and imaginary-time correlators, and understand the origin of the imaginary 
part found in [23-27]. 

The discussion will be general and pedagogical: our purpose is to illustrate 
simple and general physical effects, rather than to obtain quantitative results. 
We shall focus on a specific correlator, which is defined in Sec. 2, where we 
also recall basic properties of real and imaginary time correlators and their 
relation through analytic continuation. Then we discuss several models for 
the medium in which the heavy quarks propagate. The medium is assumed 
to be a non relativistic cold Fermi gas in Sec. 3, a thermal bath of mesons 
in Sec. 4 and a thermal bath of photons and charged particles in Sec. 5. 
The interaction between the two heavy fermions is an instantaneous potential 
in the first case, is mediated by meson exchange in the second case, and 
photon exchange in the latter. We shall see that in the limit where the mass 
of the heavy fermions becomes infinite, the correlator obeys indeed, for large 
times, a Schrodinger equation with an effective potential that accounts for 
infinite resummation of static interactions. In Sec. 5, we show that the effective 
potential can develop an imaginary part corresponding to scattering processes 
between the heavy quarks and the light fermions of the thermal bath. Finally, 
in Sec. 6 we summarize our results. 



2 Some definitions 

We shall consider a pair of heavy fermions, denoted respectively by Q and 
Q, which are supposed to be sufficiently heavy that they can be described 
by non relativistic quantum mechanics. We shall denote the quantum fields 
associated to Q and Q respectively by ip and x- We shall often refer to Q and 
Q as to a quark and an antiquark, although, in the simplified models that we 
shall discuss, these will not necessarily have all the attributes of the quarks of 
QCD. 

We shall focus our discussion on the evaluation of the following propagator: 

G > (t,r 1 ;t,r 2 \0,r[;0,r / 2 ) = ( X (t, r 2 ty(t, r^O, r'Jx^O, r' 2 )), (1) 

which represents the probability amplitude to find, at time t, a pair of heavy 
quarks at positions r , 1 ,r 2 once such a pair had been added to the system at 
time t — at positions r[, r' 2 . We also define 

J Q (t; r u r 2 ) = r 2 )i/>(t, n), f Q (t; r u r 2 ) = ^(t, rJxHt, r 2 ), (2) 
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with the time dependence given by the Heisenberg representation 

J Q (t; ri ,r 2 ) =e im J Q ( ri ,r 2 ) e~ iHt . (3) 

Here, H is the full hamiltonian, which can generally be decomposed into three 
contributions: 

H = H Q + H med + H int , (4) 

where Hq is the (non relativistic) hamiltonian describing the heavy fermions 
in vacuum, H me d is the hamiltonian of the medium in which the QQ system 
propagates, and H int represents the interactions between the medium and the 
heavy fermions. Various model hamiltonians will be considered in this paper. 

In general, the expectation value in Eq. (1) will be a thermal average, 

G>(t, r i; t, r 2 |0, r\- 0, r' 2 ) = ^Tr {e^ H J Q (t; r u r 2 )jj(0; r[, r' 2 )} , (5) 

with Z = Tre~ /3// and (3 the inverse temperature, (3 = 1/T. If we denote by \n) 
the eigenstates of H, and by E n the corresponding eigenvalues, one can write 
for G > (t) the following formal 1 expansion (in order to alleviate the notation, 
we omit the spatial coordinates): 

G > (t) = jJ2 e ~ PEn T, ei{En ~ Em)t ( n \ J Q\ m )( m \ J h\ n )- ( 6 ) 

n m 

As suggested by Eq. (6), G > (t) is an analytic function of the (complex) time 
t in the strip —f3 < Im t < 0. For t = —i(3, this function takes the value 

G>(t = - i( 3) = ±Tr{j Q e-^4}. (7) 

We shall return to this expression shortly. 

One can introduce, together with G > (t), a collection of related correlators (see 
e.g. [28]; the conventions employed in this paper are those in [29]). Thus, one 
defines 

G<(t, n; t, r 2 |0, ri; 0, r' 2 ) = (4(0; r[, r' 2 )J Q (t; n, r 2 )>, (8) 
related to G > (t) by the KMS condition: 

G < (t) = G > (t-ij3), (9) 

and the retarded propagator 

G R (t)=i9(t)[G>(t)-G<(t)}, (10) 
1 Formal, because we treat all states here as discrete states. 
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as well as the time-ordered propagator 

G(t) =ie(t)G > (t)+ie(-t)G < (t), (11) 
which is used in perturbation theory. Note that 

G(t) = G R (t)+iG < (t). (12) 

From the Fourier transforms G > (uj) and G < (uj) one obtains the spectral func- 
tion p G : 



/oo 
dte iuJt [G>(t) - G < (t)} = (1 - e-^)G > {uo) = (e^ - 1)G < 
-oo 

(13) 

the last equalities following from the KMS condition G > (uj) = d 3uj G < (uj). One 
can use the spectral function to define an analytic propagator (z complex): 

G(z) is an analytic function of z in the whole complex plane except the real 
axis. When z approaches the real axis from above, G goes over to the Fourier 
transform of the retarded propagator: G r {uj) = G(u + irj), with uj real. The 
analytic propagator is also related to the propagator that enters perturbative 
calculations at finite temperature in the imaginary time formalism, namely the 
Matsubara propagator: G(iuj n ) where uj n = 2nii/f3, n integer, is a Matsubara 
frequency. 

As to the Fourier transform of the real-time propagator, it reads (see Eqs. (12) 
and (13)): 

n( \ f +o ° d( i PG(q°,q) , . , x n ,x 

G{u,q)= ~ 5 - + ip G (u,q)N(u), 15 

J-oo 2tt q — uj — it] 

with 

Ar M = ^n- (is) 

It is useful to observe that the limits uj — > are in general distinct for the 
anaytic and the real time propagators. Using the fact that the spectral density 
Pg{%) vanishes when q Q = 0, one can write for the real time propagator 

ni n ^ f +ao dq° p G {q°,q) Pai^q) 

G(u = 0,q)= + tThm , 17 

while only the first term in the equation above is present in G{oj = 0), as is 
obvious from Eq. (14). 
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The relations above are valid not only for the propagator G, but also for most 
of the propagators that we shall encounter in this paper, to within signs for 
fermions. Thus for instance, the heavy fermion propagator reads (dropping 
again the spatial coordinates to alleviate the notation) 

s(t) = i(T^{t)^{o)) = ie(t)s > (t) -ie(-t)s < (t), (is) 

which can be also expressed as: 

S(t) = S R (t)-iS < (t) (19) 
with S R (t) = i6(t) [S>(t) + S<(t)}. 

We shall, most of the time in this paper, consider the case where the mass of 
the heavy quarks is taken to infinity. In this limit two important simplifications 
occur. First, the statistical factors describing the probability to find a heavy 
quark in the heat bath vanish exponentially (as e~ M ^ T ). It follows that the 
states involved in the thermal averages may be considered as vacuum states 
for the operators ip and \i so that S < and G < vanish. In other words, in this 
limit, the time ordered propagator and the retarded propagator coincide. The 
same situation may be imposed in the case where M is finite, by treating 
the heavy quarks as "test particles", that is as particles which are not part 
of the heat bath. The second simplification which occurs in the infinite mass 
limit is that the heavy particles do not move, so that their propagators are 
proportional to delta functions of the coordinates. To see that, recall that, 
as we just emphasized, for a test particle, or in the limit of a very massive 
fermion, the time ordered propagator is identical to the retarded propagator. 
In the absence of interactions one has then S°(t) = S^t) with 



Q0(t T )-f — f dk P -™t ik- x ~ l -jO(A f dk r -i^jt+ik- X 
S *V> X >-J 2 J(27r)3 6 6 u-jL + irj-WlJ (2^ 

1 /2M7T\i ,W» 



=m W?\~) e ' T ' (20) 

In the limit M — >oo, 

S° R (t,x)=i6(t)6(x), (21) 
so that (with G° the propagator (11) in the absence of interactions) 

lim G°(t, r i; t, r 2 |0, r[; 0, r' 2 ) = i9(t)S( ri - r[)5(r, - r' 2 ). (22) 

As we shall see, such delta functions can be factorized also in the presence of 
interactions. 

The physical information that can be extracted from the correlator (1) de- 
pends on whether one follows the evolution in real time or in imaginary time 
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(Euclidean formalism). We shall consider only the large M limit here. In 
real time, it is convenient to consider the correlator (1), or rather the re- 
tarded one in Eq. (10), from the point of view of linear response. We con- 
sider the medium initially in thermal equilibrium at temperature T, and 
perturb it at time t — with a perturbation proportional to the operator 
•/q^'i, r*' 2 ) + JQ(ri,r' 2 ). That is, for times t > 0, the hamiltionian of the full 
system is H — \[J\(r[, r' 2 ) + Jq(t ,/ 1 , r' 2 )]. The expectation value of the operator 
Jq(vi, r 2 ) + Jq(vi, r 2 ) at time t >0 is then given by the retarded propagator 

( J f (r-i, r 2 ) + J(n, r 2 )) t = A f dt'Gnit; r u r 2 \t>; r[, r' 2 ). (23) 

Jo 

In imaginary time, lattice calculations provide information on the correlator 

G>(r, r i; t, r 2 |0, ri; 0, r' 2 ) = ( J(r; n, r 2 ) J+(0; ri, r 2 )> (24) 

from which one can, in principle, reconstruct the spectral function, and then 
real time information. When the imaginary time r is set equal to f3, the phys- 
ical information that one can extract concerns essentially equilibrium proper- 
ties. Consider Eq. (7) in the limit M — > oo. It can be written as 

G>(t = -i(3) = y^(n\J Q e^ H J^ Q \n), (25) 

where H = H me d + Hint but, according to the remarks made earlier, the states 
\n) can be chosen as eigenstates of the medium Hamiltonian H me d (Hq\ti) = 
Jq\u) = 0). The sum over the states \n) in Eq. (25) can then be interpreted 
as the partition function for a system of two heavy quarks at fixed positions 
in a medium at temperature T; a basis of states for this system is indeed 
provided by states of the form Jn\n), the operator jQ(ri,r 2 ) creating a pair 
of heavy quarks at positions r , 1 ,r 2 , and |n) being an eigenstate of H med . Up 
to a product of delta functions that takes care of the normalization of such 
states, {n\ Jq{t\, r 2 ) Jq^, r' 2 )\n) = S(ri — r' 1 )S(r 2 — r' 2 ), this partition function 
is related to the free energy of the heavy quark pair in the medium. Specifically, 

G>(t = -t(3) = 5( ri - r;)5(r 2 - r' 2 ) e - pAF QQ {ri ~ r2 \ (26) 

where AFqq(vi — r 2 ) is the free energy of the heavy fermions located at 
positions r*i and r 2 , relative to the free energy of the medium in the absence 
of the heavy fermions, the latter being given by —TlnZ. This formula is 
analogous to that established in [2] which gives the change in the free-energy 
resulting from the addition of a heavy quark-antiquark pair in terms of the 
correlator of the associated Polyakov lines. 

The correlator (1) studied in this paper is related to correlators that are di- 
rectly involved in the determination of observables. In particular, it is closely 
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related to the following correlator 



G^(t,x) = (J M (t,x)Jl(0,O)), 



(27) 



where the operator J M (t,x) = q(t, x)F M q(t, x) excite mesons with quantum 
numbers that depend on Tm- In particular, in the vector channel, such a 
correlator enters the calculation of the dilepton production rate, a quantity 
of experimental interest, used for instance as a diagnosis of the the state of 
matter produced in the ultra-relativistic heavy-ion collisions. The correlator in 
Eq. (27) is also the quantity which is evaluated for imaginary times in lattice 
simulations, and used to reconstruct the mesonic spectral functions. Finally, 
a correlator such as (1) can also be used as a starting point for a T-matrix 
approach such as that developed in Ref. [3]. 

The correlator (1) is not a gauge-invariant quantity, since the QQ pair is 
not produced by a local current operator. Although this is in principle not 
a problem it may be preferable, when one performs approximations, to work 
with gauge invariant objects. This may be achieved by introducing a modified 
propagator in which a parallel transporter is introduced to connect the heavy 
fermions at times t = and t, as was done for instance in Ref. [23]. However 
for the purposes of the present paper, namely the study of the QQ free-energy 
and the large (real) time behavior of the heavy fermion correlator, this is not 
needed, as we shall verify explicitly in Sec. 5. 



3 Non relativistic setting at T = 

We start with the simple case where the two heavy fermions propagate in 
vacuum and interact through a spin- independent instantaneous potential V. 
Although we are dealing here with a problem of non relativistic quantum me- 
chanics, we shall use a field theory formulation that will make forthcoming 
generalizations easier. Thus we write the hamiltonian Hq in second quantiza- 
tion (we set H— 1): 



G > (*,r 1 ;t,r 2 |0,ri;0,r' 2 ) = (0| X (t, r 2 )^(t, n)^(o, ^(0, r' 2 )|0>, (29) 



where |0) is the vacuum of the operators ip and x ("0 1 0) = x|0) = = 




(28) 



and the correlator in Eq. (29) is 



(01^ = 0). 
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From the equation of motion of the fields in the Heisenberg picture 

idtij>(t, n) = T(n)^(t, n) + J dy x j (t, v)v( ri - y) X (t, y)^(t, n) , (30) 
and similarly for X , one deduces the evolution equation for G > : 

(id t - T( n ) - T(r 2 )) G>(*, n; f, r 2 |0, r' i; 0, r' 2 ) 

= 1^ (0|^(f,*)y(*-r 2 Mf,aOx(t, r 2 )V(t,r 1 )V t (0,r / 1 ) X t (0,r 2 )|0) 

+ / cfy <0| X (t, r 2 )x f (i, y)V(rx - y)x(t, vMt, ^(0, r'Jx^O, r 2 )|0), 

(31) 

where T(r) = — V^/2M. In general, this would not be a closed equation, its 
RHS involving a six-point function. However, the first term in the RHS of 
Eq. (31) does not contribute since (0\^(t,x) = 0; the last one can be easily 
evaluated by employing the equal-time anti-commutation relations 

{x(t,r 2 ),x%y)} = S(r 2 -y), (32) 

and it can be expressed entirely in terms of G > . One is left then with the 
closed equation 

(id t ~ T(n) - T(r 2 ) - V(n - r 2 )) G > (t, n; t, r 2 |0, r[; 0, r' 2 ) = (33) 

which, not unexpectedly, has the structure of a two-particle Schrodinger equa- 
tion with the boundary condition G > (0, r±; 0, r 2 |0, r[; 0, r' 2 ) —8(ri— T' 1 )5(r 2 — 
r 2 ). In fact, denoting by ^/ m (ri,r 2 ) the stationary states of the Schrodinger 
equation for the two heavy particles interacting with the potential V, we have 

G>(t, n; t, r 2 |0, r' i; 0, r' 2 ) = £ e- i£ -*(0|x(r- 2 )V(r 1 )|m)(m^ t (r / 1 )x t (r- 2 )|0) 

m 

= Ee _i ^ t *m(ri,r 2 )^(r / i y 2 ). (34) 

m 

It is useful, in view of the forthcoming discussions, to write the solution of 
Eq. (33) as an expansion in powers of V. In order to do so we note that the 
time-ordered propagator is simply G(t) = i6(t)G > (t) since G < (t) = here. 
One can then get the desired expansion for G > (t) from the usual perturbative 
expansion for the time-ordered propagator: 

00 (—j) n r+oo r+oo 

- iG(t, n; t, r 2 |0, ri; 0, r' 2 ) = £ LJ_ / ^ / 

n _Q J -co J -co 

x (OlT^^.-.^^x^r^a, ^(0,^(0,^)) |0>, (35) 
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Fig. 1. The ladder diagrams contributing to the two-particle propagator. The wavy 
line denotes the instantaneous interaction potential. 

where the fields are in the interaction picture and the interaction hamiltonian 
reads: 

H^x ) = e iH ^ Hie -iH x = JdxJ d A y^(x) X \y)U(x - y)x(y)i/>(x) , (36) 
with: 

U(x -y) = V(x - y)5(x° - y°) , (37) 

and Hq is the heavy quark kinetic energy. Note that, since (0\H = H\0) = 
0, the various contributions to the series in Eq. (35) vanish unless all time 
arguments x° are between and t. These contributions correspond to the 
ladder diagrams depicted in Fig. 1. The first term in the expansion is simply 
the product of two non-interacting single-particle propagators 

-iG°(t, n; t, r 2 |0, r\- 0, r' 2 ) = -S°(t, r x - r[) S°(t, r 2 - r' 2 ) . (38) 

In the limit M — > oo, it is easy to resum the series of ladder diagrams in 
Eq. (35) in a compact form: the intermediate propagators are proportional 
to delta functions of the positions of the particles (see Eq. (21)), and the 
integrations over the time variables can then be trivially done. The resulting 
series reduces to that of an exponential. One gets: 

G(t, n; t, r 2 |0, ri; 0, r' 2 ) = G°(t, n; t, r 2 |0, r[; 0, r' 2 ) G(t, n - r 2 ), (39) 

where G° is given by Eq. (22) and 

G(t, ri -r 2 ) = exp[-i^(r 1 -r 2 )] . (40) 

(Note that G is essentially G>: G > (t,r 1 ;t,r 2 \0,r' 1 ;0,r' 2 ) = 5(r l - r;)5(r 2 - 
r' 2 )G(t, ri—r 2 ).) This result can of course be obtained directly from Eq. (33), 
noting that in the infinite mass limit, the kinetic energy vanishes so that the 
coordinates ri, r 2 play the role of parameters in the differential equation. 

One can also perform easily the analytic continuation to imaginary time and 
get the Euclidean correlator. In particular, setting it — (5 in Eq. (40), one 
obtains that the change in the free energy is independent of the temperature 



10 



and equal to the potential V(r\ — r 2 ) (see Eq. (26)). Note that the same ana- 
lytic continuation performed before taking the infinite mass limit (i.e., setting 
* = -%P in Eq. (34)) yields 

G >(t = -ifi, ri ; -iP, r 2 \0, r[; 0, r' 2 ) = £ e^^n, r 2 )^(r' 1? r' 2 ). (41) 

In particular, for r x = r\ and r 2 = r' 2 , 

G>(-iP, n; -2/5, r 2 |0, ri; 0, r 2 ) = ]T e -^|* m (n, r 2 )| 2 (42) 

is proportional to the probability to find the two heavy particles at positions 
r*i, r 2 when these are in thermal equilibrium at temperature 1/p. 

We turn now to the case where the heavy particles propagate in a medium 
that they can polarize. We model this medium by a set of of light fermions 
of mass m, described by the field 0. These light fermions interact with the 
heavy particles with the interaction potential V. The hamiltonian takes then 
the form in Eq. (4) with Hq defined in Eq. (28), 

H med = J dx(f>\x) 4>( x ) - \ jdxdy f {x)<j>\y) V(x - y)(f)(y)(f)(x), 

(43) 

and 

H^t = - Jdxdy tp\x)(f)\y)V(x - y)(f)(y)^(x) 

+ Jdxdyx ] {x)(t>\y)V{x-y)4>{y)x{x). (44) 

The correlator of the heavy fermions is still given by Eq. (29) with now |0) 
the ground state of the many body system, which however still plays the role 
of the vacuum for the heavy fermions. 

The exact evolution equation for G > reads now: 

(id t -T x -T 2 -V( ri -r 2 )) G>(t, r i; t, r 2 |0, r[; 0, r' 2 ) = 

J dx{0\<j>% x)V(x - r 2 )0(t, x) X (t, r 2 )4>(t, n)^t(0, r'Jx^O, r' 2 )|0) 

- J dy(0\ X (t, r 2 )0t( t , y) V ( ri - y)</>(t, y),f>(t, n)^(0, r'Jx^O, r' 2 )|0) . (45) 

This equation reveals the full complexity of the many-body problem that we 
are dealing with. It cannot be solved exactly, nor in general be reduced to a 
simple Schrodinger equation, or a closed equation for G > . However the struc- 
ture of the corresponding time-ordered propagator can be easily analyzed in 
perturbation theory, as we did before. Typical diagrams are displayed in Fig. 2, 
where several types of processes are exhibited: self-energy corrections, and var- 
ious modifications of the effective interaction between the heavy fermions. In 
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fe) (b) (C) 

Fig. 2. Some of the processes occuring when a of a QQ pair propagates in a medium 
of light fermions. The diagrams (a) and (b) are accounted for by the effective inter- 
action. On the contrary diagrams (c) is not included in the effective interaction. 




VWAAj = 'WWV + 

Fig. 3. The diagrams contributing to the effective interaction. 

order to progress in our analysis, we shall assume, as done in many studies, 
that the dominant effect of the medium is to modify the interaction poten- 
tial through screening effects. We can then write a resummed perturbative 
expansion by employing the following effective interaction hamiltonian (in the 
interaction picture): 

Hf(x ) = [dx [ d 4 y^(x)x\y)U QQ (x-y) X (y)^) + 



\J dxj d 4 y^(x)^(y)U QQ (x-y)^(y)^(x) + 

^jdxj d 4 y X ^x) X ^y)U^( x - y)x(y)x( x ) ■ (46) 



+ 

2 J 

+ 



The diagrams contributing to the effective QQ and QQ (and Q Q) interactions 
are displayed in Fig. 3. The effective interactions are conveniently expressed 
in terms of their Fourier transfoms: 



U (uj,q) 



V(q) 



1 - U(u, q)U (u, q) 1 - U(co, q)V(q) 



Uo(u,q)IL(u,q)U (u,q) 
1 - U(u,q)U {uJ,q) 



(47) 
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Fig. 4. An example of the diagrams included in the exact exponentiation occurring 
in the M^oo limit. The wavy lines with black dots denote the effective interaction 
(see Fig. 3). 

where U(ui, q) is the particle-hole bubble. Note that the independence of IT 
makes the medium induced interaction no longer instantaneous. 

We can now study how the medium-modified interaction affects the two- 
particle correlator, though being aware that this does not exhaust all the 
possible processes occurring during the in-medium propagation. Note that in 
the case of a long-range interaction, like in the electron-gas problem, this choice 
is sufficient to account for most medium effects [30] . However, our purpose here 
being only to illustrate some specific effects, no further effort will be made to 
justify keeping only this particular class of diagrams. In the same spirit we 
shall ignore tadpole diagrams, assuming that the medium has the necessary 
properties to ensure their cancellation. A typical diagram contributing to the 
two-particle propagator is then given in Fig. 4. 

The treatment of the effective interaction remains a difficult task in general, 
but can be done exactly in the limit M — > oo. One starts from the pertur- 
bative expansion given in Eq. (35) with the effective interaction hamiltonian 
introduced in Eq. (46). In the M — > oo limit the various contractions that 
involve a given heavy fermion contibute an overall factor 0(t)5(r — r'). Hence, 
in any order, one can factorize 6 , (t)5(r 1 — r[)S(r 2 — r' 2 ) corresponding to the 
non-interacting result given in Eq. (22). The different ways of contracting the 
fields simply amount to the different temporal orderings of the interaction ver- 
tices x® and The integrations over these variables can be done freely and 
lead again to an exponentiation. The final result takes the form of Eq. (39) 
with here 



G(t,r 1 -r 2 ) = exp 



-if dx°f dy°{U QQ (x -y°,r 1 -r 2 )+U QQ (x -y ,0 

(48) 

By expressing the effective interaction in terms of its Fourier transform, one 
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obtains: 



G(t,ri - r 2 )=exp 



r duj r dq 1 — cos(c<jt) 
2 7 2W J2nf X 



x 



(U qq (uj, q) + e iq < ri ~ r ^U QQ {uj, q 
The large time behavior of this expression is easily obtained by using 



(49) 



One gets: 
where 



lim l-cosM 



G{t,r 1 -r 2 ) ~ exp (-itV 00 (r 1 - r 2 )) , (51) 



Ko(r-i -r 2 )=J-^ (u QQ (cu = 0,q)+e^-^U QQ (u = 0,q)) . (52) 

Quite remarkably, aside from the self-energy corrections, this large time be- 
havior (51) is the same result as would have been obtained with an instanta- 
neous potential Kx>(^i — r 2 ) = U QQ (oo = 0,r 1 — r 2 ) (compare with Eq. (40) 
that expresses the resummation of ladder diagrams involving the instanta- 
neous potential V(r\ — ^2))- We are facing here an issue closely related to the 
problem of recovering a potential model from a two-body relativistic equation 
(with the energy dependence of the interaction arising from the propagator 
of the exchanged meson). Also in this case, in order to get, in the large mass 
limit, an energy- independent interaction kernel, one has to sum all possible 
ladder and crossed diagrams [31]. Note that the evolution equation obeyed by 
G{t) at large time is simply: 



lim 



id t -J (u QQ (u = 0,q) +e l i< r ^U QQ (uj = 0, q 



This is the analog of Eq. (4.2) in [23]. 



G(t, ri -r 2 )=0. 

(53) 



So far, we have been working in real time at T = 0, considering the propagation 
of the QQ pair in a cold-dense medium of light particles. However, one can 
easily extend the previous calculations to address the case where the heavy pair 
propagates in a medium at finite temperature. This can be done by exploiting 
the analyticity properties recalled in Sec. 2, or by repeating the calculation 
above in the imaginary time formalism. In fact, when it is carried in the 
imaginary time formalism, the preceding perturbative analysis goes through 
unchanged except for the replacements t — > —it, —iU(t = —it) — > U{t) 
(in Fourier space, the relation between U(cu) and U{uj) is the same as that 
between, respectively, the real time and the analytic propagators discussed in 
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Sec. 2). In the infinite mass case, the exponentiation of the interaction effects 
still holds, and one gets 



G{-ir,r 1 -r 2 ) 



= exp 



dr'j^ dr" (U QQ (t' - r", r i - r 2 ) +U QQ (r' - r ", n - r 2 



• (54) 



At this point, it is convenient to express the interaction in Fourier space, for 
instance: 

Z^V-T>i-r 2 ) = ^ e "^ V) / 7^%^ (ri " r2) W Q<3 (^ n ,q) (55) 

P n 7 ( 27r ) 



where u; n = 2nirT is a Matsubara frequency. The integals over r' and r" in 
Eq. (54) can then be easily done. In particular, for r = (3, one gets 



G{—if3, ri-r 2 ) = exp 



-0 



(2^ 



D iq-(r 1 -r 2 ) 



(U QQ {uo = 0, <?) + W QQ (cu = 0, <? 



allowing us to identify the free-energy as 



(56) 



I 



dq 



(U QQ {u = 0, q) + U qq (uj = 0, qfj . (57) 



Note that, in spite of the similarity between Eqs. (57) and (52), there is a 
difference between the two expressions, of the same nature as the difference 
between the static limits of the real time and the analytic propagator that has 
been already emphasized in Sec. 2. Eq. (57) corresponds to the static limit of 
the analytic propagator and is always real, as appropriate for a free energy. In 
contrast, V^, given by Eq. (52), may contain an imaginary part. An explicit 
example will be provided in Sec. 5. 



We have addressed in this section some limiting cases in which it was possi- 
ble to get a closed, Schrodinger-like, evolution equation for the two-particle 
correlator: the trivial case of in-vacuum propagation and instantaneous inter- 
action; the case of a retarded in-medium interaction between infinitely-heavy 
fermions. In the latter case, the result follows from an exponentiation property 
that holds in the limit of a large mass for the heavy fermions, and asymptot- 
ically for large times. The following section will clarify the origin of the expo- 
nentiation by using an explicit model for the interaction, based on a simple 
meson exchange. 
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to- iP 

Fig. 5. The complex-time contour C = C + U C_ U Co used in the calculation of the 
real time QQ correlator at finite temperature. 

4 Interactions due to meson-exchange 



We now consider a model in which the heavy fermions propagate in a thermal 
bath of scalar mesons of mass m, described by the field $, to which they are 
linearly coupled. The corresponding hamiltonian reads: 



H = J dx^(x)(^ + g^(x))ij(x)+l dx x t (a .)|_X_ + g ^ x ) \ x ( x )+ 



-V 2 



+ / dx 



\n\x) + \{V^x)f + l -m^{x) 



, (58) 



where ir(x) is the conjugate momentum of the field <&(x). This hamiltonian 
is of the same form as in Eq. (4), with Hq the kinetic energy of the heavy 
fermions, H mei i the hamiltonian of free scalar mesons, and H int the linear 
coupling between the fermions and the scalar field. Note, however, that the 
dynamics of the two heavy fermions in vacuum, involves the full hamitlonian: 
in vacuum, as well as at finite temperature, the heavy fermions interact with 
each other through meson exchange. 

We want to evaluate the correlator 

G>(t, n; t, r 2 |0, r' i; 0, r' 2 ) = ^Tr {e^ H x (t, r 2 )^(t, n)^(0, ri)x f (0, r' 2 )} , 

(59) 

with Z = Tre~ pH . This can be done either by working in imaginary time, 
and then perform an analytical continuation to real values of t by exploiting 
the analyticity properties of G > (t) or by using Keldish-Schwinger countour 
techniques, such as the one depicted in Fig. 5. We shall use the latter in what 
follows, returning to the analytic continuation towards the end of this section. 
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The formalism used in this section, and the following one, follows closely that 
employed in Ref. [32] for the discussion of the damping rate of a hard fermion. 
Our strategy to calculate G > will rely on the following expression for the 
(time-ordered) propagator: 

G(t, n; t, r 2 |0, r[; 0, r' 2 ) = Z' 1 J [P$]G*(t, n; t, r 2 |0, r' l5 0, r' 2 ) (60) 

where G$ is the pair propagator in a given background field &(x,t) and 

Z = J[D<S>]e iSc[ * ] . (61) 

The propagator G^(t,r 1 ;t,r 2 \0,r' 1 ;0,r' 2 ) is simply the product of the two 
propagators of the heavy fermions in the external field $: 

iG*(f , rn f, r 2 |0, r' i; 0, r' 2 ) = S*(t, n - ri) S*(f , r 2 - r' 2 ). (62) 

The average over the configurations of <£>, repesented by the functional integral 
in Eq. (60), can be done easily since the action is quadratic: 

S c m = ~lJ c d 4 x J c d*y*(x)Dc\x-y)*(y). (63) 

The field $ is defined on the contour C in the complex-time plane displayed in 
Fig. 5 and obeys the periodicity condition $(£) = Q(t — if3). The corresponding 
propagator along the C contour satisfies the KMS conditions and is given by: 

-iD c (t) = d c (t)D>(t) + e c (-t)D<(t), (64) 

In the above, in order to give an unambiguous meaning to 0c(t), one implicitly 
assumes the path in the complex-time plane being parametrized by a real 
parameter u, t = t(u), with u increasing as t runs along the contour. Note 
that the propagator along C + coincides with the real-time propagator defined 
in Eq. (11), while the one along Cq is closely related to the imaginary-time 
Matsubara propagator, namely 

-iD Co (t = -tr) = A(r). (65) 

The usefulness of Eq. (60) relies on the possibility to calculate explicitly 
which can be done easily in the limit M — > oo, to which we shall restrict 
ourselves here. In this limit the propagator S$ coincides with the retarded 
propagator (note that in Eq. (60) S$, or G$, satisfies the same boundary con- 
ditions as G, namely retarded conditions; these differ from the KMS conditions 
satisfied by the meson propagator). We have: 

S*(t, r t - r-0 = i0(t)6( ri - rj) exp (-ig J* dt'$(t', r<)) , (66) 
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so that 



G*(t, r i; f , r 2 |0, r' l5 0, r' 2 ) = i^^n-r^r,-^)^, n - r 2 ), (67) 
where 



and 



G<s,(t,r 1 - r 2 ) = exp (-i J d 4 zJ(z)$(z)^j , 

J(z) = g f dt'5(z° - [S(z - ri ) + S(z- r 2 )] 
Jo 



(68) 



(69) 



The average over the configurations of the field $ can now be performed, using 
Eq. (60). One gets 



G(t, n - r 2 ) = exp 



f c d A x j c d 4 yJ(x)D(x-y)J(y) 



(70) 



where only the C + part of the contour contributes, and accordingly D is the 
real-time propagator of the scalar meson. Inserting the current (69) into Eq. 
(70) one gets: 



ig 2 ds ds' (D(s-s',0)+D(s-s',r 1 -r 2 )) , (71) 



G(t, r 1 - r 2 ) = exp 



which is analogous to Eq. (48) obtained in the previous section. Everything 
follows then in close analogy with what was found in the non-relativistic model 
of Sec. 3, with the meson propagator playing now the role of the effective 
interaction. In particular the large-time behavior of the QQ propagator is 
given by Eq. (51) in terms of the static propagator which, at his stage, is 
independent of the temperature: 



V^ri ~ r 2 ) = -9 2 jj^y ^ + e^~^) D{u = 0, q) 

= ~9 2 [ 7^% (l + e^~^) — L^. (72) 
J (27r) d v ' q 2 + m 2 

The effective potential is attractive, as expected for scalar exchange. The first 
term, which represent a self-energy correction, can be evaluated in dimensional 
regulatization (as done e.g. in Ref. [23]), while the second one yields a Yukawa 
potential. One gets (with r=\ri — r 2 |) 

V„(r 1 -r 2 ) = £L-tpl. (73) 



We now allow the scalar meson to couple with massless fermions through the 
same Yukawa coupling. The thermal bath is then composed of scalar mesons 
described by the field $ and massless fermions. On the other hand the heavy 
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fermions are still treated as test particles. Under the hypothesis that these 
heavy fermions interact mainly with the soft (i.e., carrying momenta qo, \q\ <C 
T) modes of the scalar field, the effective action for the averaging over the $- 
configurations will be given by the "hard thermal loop" (HTL) effective action 
of the soft modes (see the corresponding discussion in Ref. [32]). This is simply 
given by the same gausiann expression as in Eq. (63), but with Dq embodying 
now also self-energy corrections accounting for the effect of the light fermions 
on the meson propagation. There is of course a close analogy with the particle- 
hole bubble summation of the previous section. In the HTL approximation the 
meson self-energy (see Appendix A) is momentum-independent and provides 
a mere correction to the mass of the meson, which becomes 



m 2 D (T) = m 2 + IT = m 2 + iV dof 



g^r 2 

6 



(74) 



with the factor A^of counting the internal degrees of freedom of the light 
fermions other than spin (flavor, color...). This increase of the mass corre- 
sponds, in coordinate space, to the screening of the interaction potential and 
m D is commonly referred to as the Debye screening mass. 

We turn now to the Euclidean propagator, where the heavy-quark current is 
defined along the C part of the contour. Eq. (71) becomes then 



G(-ir,r 1 -r 2 ) = exp 



(2k) 



(75) 

where the imaginary-time propagator is defined in Eq. (65). By expressing the 
latter in Fourier space 



A(t'-t", q) = - £ e-^'-T") A{iuJn , q) > 

P n 

one can easily performs the imaginary time integrals and get 

d( l , iq-(ri-r 2 ) x 



(76) 



G(-ir,ri - r 2 ) 



exp 

.2 



g 



f—( 

J (2tt) 3 V 



l+e l 







x 



— A(iuj n = 0,q) + - 2(1 - cosu n T)A(iu n , q) 

where u n = 2nirT is a Matsubara frequency. In particular, for r 
cosine does not contribute and one has simply 



G(-i(3,r 1 - r 2 ) = exp 



fig' 



dq° p D ( 



(2tt 



2tt 



, (77) 
P, the 

, (78) 



where we have used the expression (14) in order to express the Matsubara 
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propagator A(iu n = 0) in terms of the spectral function pu of the meson 
propagator. 

Note that although the expression (78) looks like the analytic continuation of 
the expression (51) (with the corresponding V^), much information contained 
in the Euclidean propagator (80) is lost when we specify r = f3, because of the 
cancellation of the terms involving non vanishing Matsubara frequencies (the 
cosine term in Eq. (77)). In the expression (78), the frequency integral is the 
static limit of the (analytic) propagator, which coincides here with the static 
limit of the real time propagator (this is not so in general, as the example in 
the next section will show). 

Further insight on this issue can be gained by proceeding slighlty differently 
so as to make the analytic continuation to real time straightforward. Let us 
then express the Matsubara propagator in Eq. (77) in terms of the spectral 
function po as follows 

A(r' - r", q) = j ^V^'^'W<?o, q)W ~ O + (79) 
This also allows us to perform the integrations over r' and r", with the result 



/ 2 f dc t ft , iq.(ri-n)\ [ d( i pD(q°,q) \ 

G(-.r, ri - r 2 ) = exp | 5 J — (1+e »>) J ^^^r| x 



dq°p D (q°,q). _ o 



(e-^-l)(l + iV(g )) x 



x 



2vr (g ) 2 j 



The same expression, with r replaced by it, holds for G(t, r± — r 2 ), as expected 
from the analyticity of G(—ir, r 1 — r 2 ) for r < (3. One can then rearrange the 
terms in order to recover the results obtained earlier, e.g. Eq. (72). However, 
in the limit t — /3 the last two exponentials in Eq. (80) cancel each other and 
one is left with Eq. (78). 

It follows from Eq. (78) that the free energy of the QQ pair is given by 

A W )^ W ^-fl£^:. ( 81 ) 

Thus in this simple model the free-energy that is obtained from the imaginary- 
time two-particle propagator at T — /3, can be identified with the potential to 
be used to study the real-time propagation of a QQ pair. In the next section 
we shall find that extra contibutions to the potential may arise, that are not 
captured by an imaginary-time equilibrium calculation at t — /3. 
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5 Heavy fermions in a hot QED plasma 



We now discuss the case, closer to hot QCD, in which the QQ pair is placed 
in a thermal bath of photons and light fermions. We consider only the limit of 
infinite mass for the heavy fermions, and, to be specific, we choose the Coulomb 
gauge to describe the electromagnetic field (our final results are independent 
of the choice of the gauge). The hamiltonian describing the system reads: 

H = gj dxA (x) (-^{x)^{x) + xHxMxj) + H Ch + H t (82) 

with ifcb the photon hamiltonian in the Coulomb gauge, and Hf the hamilto- 
nian of the light fermions, including their interaction with the electromagnetic 
field. 

In order to evaluate the real-time propagator of Eq. (1) we proceed as in the 
previous Section and consider first the propagation of the heavy fermions in a 
given background field A^. Then we average over all the possible configurations 
of the gauge field. The treatment given here is similar to that of the fermion 
damping rate presented in Ref. [32]. One gets, using the notation of Eq. (67): 

G(t, ri - r 2 ) = (exp (-i J d 4 zJ (z)A (z)^j) , (83) 

with 

J (z) = f dt'5(z° - f) [-gS(z - n) + gS(z - r 2 )\ . (84) 
J o 

Assumimg that the heavy fermions interact mostly with the soft modes of the 
electromagnetc field, we use the HTL effective action in order to perform the 
average over the gauge-field configurations. One obtains 



G{t,r x - r 2 ) = exp 
= exp 



% -\ d A x d 4 y J (x)D 00 {x -y)J (y) 
ig 2 f ds f ds' (Doo(s - s',0) - Doq(s - s',r 1 -r 2 )) 

JO JO 

(85) 



where D^ v (x — y) is the real-time HTL photon propagator satisfying the KMS 
condition. The large-time behavior of the evolution equation is still given by 
Eq. (51), with related to the static propagator, as in Eq. (72). The real-time 
photon propagator in Fourier space is given in Eq. (B.7) and, for uj — 0, 

— 1 7im 2 n T , . 

fl ' ,u = MI = fT4 + 'sf<P' (86) 
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The effective potential reads then: 



(27T) 



7rm 



T 



qr 2 + ra|> 



\q\{q z + mjj) 



2 \2 



, (87) 



which agrees with that found in [23]. Note that, as it should, the (real part of) 
the potential between the quark and the anti-quark is attractive and screened. 
Furthermore, contrary to what happened in the scalar-meson model of the 
previous section, the potential develops an imaginary part. This imaginary 
part originates from that of the HTL photon self-energy (see App. B) and, 
physically, results from the collisions between the light fermions of the hot 
medium and the heavy quarks (see App. C and below). Apart from a difference 
in the sign of the self-energy contributions, the real part is identical to that 
found in the scalar case. The angular integration in the imaginary part can be 
easily performed and one is left with (adopting the same notation as in Ref. 
[23]): 



m D + 



i——(p(m D r) , 

47T 



where r = \r 1 — r 2 \ and 



4>{x) 



poo 

2 / dz- 



z 2 + iy 



sia(zx) 



zx 



(88) 



(89) 



Thus G behaves as a damped exponential 
G(t, ri 



r 2 J ~ exp 



g 2 ( e-^DT" 

i— m^H 

47r V r 



exp 



An 



(90) 



The damping factor increases with the temperature, as shown in Fig. 7: the 
function <f>(x) vanishes for x = and increases monotonously, approaching 1 
as x — > +oo. In the r — > oo limit the contributions of the two heavy fermions 
to the pair propagator factorize and, accordingly, the damping factor in Eq. 
(92) results from the sum of two equal terms (for Q and Q respectively): 



lim <f>(x) 



1Q = 1Q Z 



8tt 



(91) 



which coincides with the v — > limit of the heavy fermion damping factor in 
Eq. (3.6) of [33], after dropping the Casimir factor (the calculation of Ref. [33] 
being done in hot QCD). Thus, the collisional damping rate is most important 
when the heavy quarks are far apart. When they get closer together, interfer- 
ence occurs with the process by which the exchanged photon gets absorbed 
by the light particles in the heat bath. When the QQ separation vanishes, this 
interference is completely destructive and kills the imaginary part. 2 



We therefore differ in our interpretation from that suggested in Ref. [27] : the "dis- 
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Fig. 6. The dependence on the QQ separation of the function 4> describing the 
damping of the in-medium two-particle propagator. 

In order to get a feeling for the magnitude of the effect of the imaginary part on 
the QQ propagation in a hot medium, we have plotted in Fig. 7 the behavior 
of the damping factor 

7 = g -^<t>{m D r) (92) 

arising from Eq. (90), as a function of the temperature. We shall consider 
the propagation of a heavy fermion pair in a hot QED plasma (identifying 
g 2 1 '4-7T = «qed) consisting of photons, electrons and positrons. In order to deal 
with some sensible numbers we take the heavy fermion pair as being separated 
by a distance equal to mean value of the radius of a ji + /i~ atom in its ground 
state, namely: 

r = (r)is = ^Bohr = « 3.89 MeV" 1 , (93) 

2 2 /iftQED 

\x being the reduced mass of the fi + fi~ pair. One sees in Fig. 7 that the damping 
is quite substantial. Actually the real growth of the damping factor with the 
temperature should be even faster, since one expects that as the temperature 
increases the effective potential gets more and more screened leading to a lower 
binding energy and a more spread wave-function. This would entail (treating 
the imaginary-part as a small perturbation) a larger value for (7)15, since 
the function in Eq. (92) increases monotonically with the QQ separation. 
A more refined estimate of the in-medium quarkonium decay-width has been 
given recently in [23]. 

Turning now to the Euclidean correlator, we note that the analytic continu- 
ation to imaginary-time can be performed by following the same steps as in 
the scalar case. One then easily obtains 

appearance" of the exchange meson is in fact acting against the major contribution 
coming from the individual damping factors of the heavy fermions. 
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Fig. 7. The damping factor for a heavy fermion-antifermion pair placed in a hot QED 
plasma (of photons, electrons and positrons) as a function of the temperature. The 
distance between the heavy fermions has been set equal to the mean radius of the 
ground state of a pT atom. 



G(-if3, ri - r 2 ) =exp 



=ex p{"* 2 /(^ (i " e " 



( ri _ r2 )\ i 



, 2 , ( 94 ) 
q 2 + mh) 



where in the last line we have employed the sum rule given in Eq. (B.8). The 
free-energy change occurring when the QQ pair is placed in the thermal bath 
is then: 

AWH-^-f^. (95) 

m D being the Debey screening mass of the photon. Note that G can also be 
written as the thermal average of two Polyakov lines (we assume the average 
to be made over Euclidean fields Af(r) such that A (— it) = iAf (r)): 

G(- t f3, ri -r 2 ) = (exv(igJ Q P dT'Af (r', ri )j exp(-tg [dr'Af (r', r 2 )j ). (96) 

We recover in Eq. (95) the expected lowest order perturbative result for the 
QQ free energy. However, comparing Eq. (95) with Eq. (90) (which refers 
to the real-time propagation problem), one can see that considering only the 
imaginary-time result at r — f3 one misses the damping contribution arising 
from the collisions with the particles of the thermal bath. 
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6 Conclusions 



We have shown in this paper that the real-time propagator of a (infinitely) 
heavy QQ pair crossing a hot medium obeys, at large times, a Schrodinger 
equation. We have evaluated the effective potential entering this Schrodinger 
equation using simple models for the medium in which the QQ propagates, as 
well as for the interaction between the medium and the heavy fermions. Our re- 
sults corroborate those of Refs. [23,24] and put them in a broader perspective. 
In agreement with [23-25], we have found that in the QED case the real-time 
potential develops an imaginary part, for which we proposed an interpretation 
in terms of the collisions of the heavy quarks with the light fermions of the 
thermal bath. We emphasized that the damping is largest when the heavy 
quarks are far apart, where it coincides with the damping of heavy fermions 
that have been calculated long ago. At small separation, interference effects 
strongly suppress the damping in the correlator. 

This effect of the collisions on the motion of the QQ, and the fate of its 
possible bound states, is distinct from the process of gluo-dissociation that 
was first considered by Bhanot and Peskin [34,35], and later by [36]. The gluo- 
dissociation reaction can be viewed as the LO pQCD result for the quarkonium 
break-up. In [11] the latter was the process considered in the study of collisional 
melting of charmonium in the QGP (namely the collision with a real gluon 
of the thermal bath). It was also shown that the above reaction could be 
considered as the analog of the deuteron photo-dissociation reaction, whose 
cross section was evaluated in [37]. The outcome of the calculation of the 
gluo-dissociation process is a width of the charmonium which decreases with 
increasing temperature: this due to the fact that, as the binding energy gets 
lower and lower, the cross section is peaked at values of the gluon energy 
corresponding to very small phase space. The origin of the imaginary part 
of the potential discussed in this paper is closer to the so-called quasifree 
destruction mechanism explored in [38] (the scattering off light partons of the 
thermal bath) and which leads to a collisional width for the quarkonium in 
the QGP that increases with the temperature [38,39]. 

A second aspect of our analysis concerns the relation between the real-time 
and the Euclidean propagators. One finds that the real part of the effective 
potental governing the large-time behavior of the heavy quark propagator can 
indeed be extracted from the Euclidean correlator calculated for imaginary 
time t = P, where (3 is the inverse temperature. In particular we found that 
the real part of the effective potential coincides with the QQ free-energy. 
However we showed that the Polyakov-line correlator is "blind" to damping 
effects giving rise to the imaginary part of the effective potential. 

A limitation of the present study is of course related to the fact that it has 
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been performed in the M — > oo limit for the heavy quarks. In [23-25] the re- 
sults obtained in this limit where directly employed to address also the finite 
mass case, plugging the potential into a Schrodinger equation that contains 
the kinetic energies of the heavy fermions. This was used in particular in 
order to reconstruct, starting from the popagator G > (t), the in-medium spec- 
tral function of the quarkonium. It should be possible however to extend the 
present analysis to treat the finite mass case using the techniques of effective 
field theories that have been developed in the T = case [40]. Note finally 
that in dealing with the M — > oo case one could drop contributions arising 
from the exchange of magnetic photons. While we expect magnetic photons 
to remain subdominant when M is finite but large, their contribution may 
generate finite temperature IR divergences (arising from the lack of screening 
for the magnetostatic modes) which can be cured by the method presented 
in [32]. 



A HTL scalar self-energy 

The scalar self-energy is defined through the Dyson equation 



written here in imaginary time, with Ui = 2irl/j3 and the free Matsubara prop- 
agator Aq 1 (iu)i, q) = oof + q 2 + m 2 . To lowest order in g, IL(iu>i, q) is given by 
the one-loop diagram involving light fermions: 



where the factor A^of counts the light fermion internal degrees of freedom other 
than spin (flavor, color). After performing the trace over the Dirac indices, 
and the sum over the Matsubara frequencies, one can perform the analytical 
continuation — > q°, with q° off the real axis. Defining p 1 =p, p 2 = p — q, 
€i=pi and e 2 =p 2 , ra(ej) = l/(e |9ei + 1), one finds: 



A-l(iui, q) = A \iuj h q) + U(iu h q) 



(A.I) 



I +°° r Art 
U( i u h q) = N doi g 2 - £ UJ^Tr[S F (iu; n ,p)S F ( 




iu n -iui,p-q)], (A.2) 




x j(eie 2 -p 1 -p 2 ) 
(e 1 e 2 -p 1 -p 2 )-^ 



(e 1 e 2 + p 1 -p 2 ) 



+ (eie 2 + p 1 ■ p 2 ) 



(A.3) 
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The Hard Thermal Loop approximation [28] consists in the following kinemat- 
ical approximation, valid when q , \ q\ < T: 



e 2 ~p-qcos9, n(e 2 ) - n(ei) ~ (e 2 - ei) 



dn 
~de~ 



~—qcos9 



dn 



(A.4) 



where 9 is the angle between p and q. In this kinematical regime the first 
and the third terms within the parenthesis in Eq. (A. 3) do not contribute. 
The remaining terms contain contributions that are UV divergent term but 
independent of the temperature or of the momentum, and can be absorbed 
into the renormalization of the meson mass. One finally gets [41] 



U(q°,q)=4N do{ g 2 J 



dp n(p) 
(2tt) 3 p 



9 2 T 2 



N, 



dof" 



(A.5) 



B The HTL real-time photon propagator 



The one-loop longitudinal photon self-energy reads (for simplicity we assume 
that the coupling with the electromagnetic field is the same for all the species 
of light fermions) : 

X +°° [dp 

n 00 (^, q) = N do{ g 2 - Yl / j7^^^[l° S F(iuJn,ph S F (iuJ n - iu h p- q)] , 

(B.l) 

leading to, after summation of the Matsubara frequencies, and analytic con- 
tinuation 



n o(<7°,g) = iW f-^h— x 

J (2ir)- s e 1 e 2 



j ,n(e 2 )-n{e l ) 1 - n(ei) - n(e 2 ) 

X 1 (ei62 + Pl " ft) g°-(ei-e a ) + (ei62 " Pl " ft) f-(e 1 + e 2 ) + 
t -» ^ n(e 2 )-n(e 1 ) 1 - n(e 1 ) - n(e 2 ) \ 

In the HTL approximation, the dominant contributions come from the " Lan- 
dau damping" process (the first and third terms in the equation above; note 
the difference with the scalar case), and we have 



n 00 (g ,qr) = 47V dof 



(2k) 



dn\ 



, i P dp 
2 'o \ de 



1 a a 1 cos6 
dcos9- 



q° — q cos 9 



~m D (l-2 ln ^T 



(B.3) 
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where x = q°/q and mjj is the Debye screening mass 

a 2 T 2 

m 2 D = 7V dof ^-. (B.4) 
Form this one deduces in particular the imaginary part (qo = w + irj) 

lmU 00 {uj + q) = -m|,^(g- (B.5) 

and, after solving the Dyson equation with the non-interacting longitudinal 
propagator Aq 1 ^ = —l/q 2 , the analytic (off the real axis) propagator 

qr 2 + m D (1 - f mf^J 

From the latter one gets the longitudinal photon spectral function which allows 
to express the real-time propagator as follows: 

Ax>(", 9) = ^ + ^ g0 _ (h; + jTy) + «Poo(", flf)JV(o;) . (B.7) 



The a; = limit of the above real time correlator can be obtained from the 
sum rule [28] 

-°° dq° Poo(q°,q) 1 i (Bg) 



d 2% q° q 2 q 2 + m 2 D 

and from the low-energy behavior: 



p 00 (tu,q)N(tu) ~ Jff 2 , 2 ~- (B.9) 
w^o IqKqr 2 + m z D ) 2 uj 



C Physical interpretation of the damping 



Here we give a physical interpretation of the processes responsible for the 
damping of the QQ correlator discussed in Sec. 5 . 

We consider a heavy fermion with four-momentum P— (M, 0) (M^>T) in a 
hot QED plasma. The interaction rate with the light fermions of the thermal 
bath is given by: 



rW = 2^/ B / B / B (*)W + fi-ft-ft)x 



x [m(l - n 2 ){\ - n 3 ) + (1 - n!)n 2 n 3 ] \M\ 2 , (C.l) 
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(C3-A) 



Fig. C.l. The collision process responsible for the imaginary part of the QQ corre- 
lator 

where f p . = J dpj (27r) 3 2ej, rii is a fermion occupation factor and the overall fac- 
tor 2 accounts for scattering over particles and antiparticles. The correspond- 
ing process is depicted in Fig. C.l: the particle 1 hits the heavy fermion and 
comes out with four- momentum P 2 . The index 3 refers to the four- momentum 
of the heavy fermion after the scattering. The second term in the parenthe- 
sis accounts for the inverse process. When the mass of the heavy fermion is 
infinite, ^3 = and the inverse process does not contribute. 

The modulus square of the amplitude, averaged over the spin of the incom- 
ing heavy fermion and summed over the spins of the other particles is easily 
evaluted [42]: 



In the above A^ is the HTL resummed photon propagator, which we take 
in the Coulomb gauge. Furthermore, in the large M limit, we can keep only 
the contribution arising from the exchange of a longitudinal photon and ignore 
that from the transverse photons. Finally, because the interaction is dominated 
by the exchange of soft momenta, e 3 ~ M and p 1 ~ p 2 - Under the above 
approximations one can write: 



The statistical factor can be rewritten as [42]: 

m(l - n 2 )(l - n 3 ) + (1 - m)n 2 n 3 = (m - n 2 )(l + N(q ) - n 3 ) . (C.4) 

In the present kinematical conditions, where p\ ~p 2 ~T, and q° = e 2 — ei <^ T, 
one can write: 



\M\* = 8g 4 A, u (q)A; x (q) [P'Pg + P^P? - g^(P-P 3 ) + g w M x 

x \P?P£ + P%P£-sr\PvP 2 )\ . (C.2) 



\M\ 2 = 8g 



A 00 (g)| 2 (2M 2 )(2e ie2 ). 



(C.3) 



(1 + N(q ) - n 3 ) 



T 



(C.5) 



,0 



and 



(m -n 2 ){l + N{q ) -n 3 ) 



dn 



(e 2 - ex) 



T 




(C.6) 



dpi 



,0 
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It is then convenient to change the integration variables as follows [42] : 




where the two Dirac-distributions account for energy conservation at the inter- 
action vertices. The first one expresses the fact that the energy transfer to the 
heavy fermion ~ q 2 /2M is negigible. The second expresses the conservation 
of energy at the vertex with the light fermions, in the HTL approximation 
where q° = e 2 — ei m V\ ■ q. One can use Eq. (C.7) to perform the angular 
integration over p l in Eq. (C.l). One is left with: 



which is precisely the contribution to the QQ damping factor arising from 
processes involving a single heavy fermion. 
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